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ABSTRACT 

The velocity law of hot star winds is usually parameterized via the so-called beta velocity law. Although this parameterization stems 
from theoretical considerations, it is not the most accurate description of the wind velocity law that follows from hydrodynamical 
calculations. We show that the velocity profile of our hydrodynamical wind models is described much better by polynomial approxi- 
mation. This approximation provides a better fit than the beta velocity law already for the same number of free parameters. 

Key words, stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics 



Hot stars lose an important part of their mass via line-driven 
winds. Consequently, mass-loss rate prescriptions are an impor- 
tant part of evolutionary codes. Such prescriptions can be de- 
rived from theoretical modelling of hot star winds. However, any 
theoretical prescription should be tested against the mass-loss 
rates derived from observations. Although there are several pos- 
sibilities for deriving wind mass-loss rates (e.g., from Ha line 
profiles, from ultraviolet P Cygni line profiles, and from radio 
continua), knowing of the wind velocity law is needed to de- 
rive the mass-loss rates from observations. Although there ex- 
ist sophisticated hydrodynamic calc ulations giving the wind ve- 
locity la w in a tabular form (e.g., iGrafener & Hamaniil l2005l 
iKrticka & Kubat 2010), analytic formulae are useful in many 
analyses. Consequently, the question of proper velocity law pa- 
rameterization is important for analysing observed hot star wind 
spectra. Here we develop a parameterization that is as close as 
possible to the results of hydrodynamical calculations. 

The most frequently used parameterization of the wind ve- 
locity law is the so-called beta velocity law. 
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where Voo, i3 is the wind terminal velocity, f3 a free parameter de- 
scribing the steepness of the velocity law, and i?* the stellar ra- 
dius. Thus, for parameterization of the wind velocity law using 
Eq. [T] we need two parameters, namely /? and the wind termi- 
nal velocity Uoo.^- The popularity of the beta velocity law Eq.[T] 
stems from its simplicity and from the belief that the velocity 
profile described by Eq. [T]is close to the real one. Also simple 
analytical wind solutions have the form of Eq. [T] For example, 
in the case of accelerating force proportional to ^ l/r^, the so- 
lution of the momentum equation (Milne 1926, Chandrase khaJ 
11934 lRublevlll965h leads to v{r) = Woo,/J v^l - im- 
ply ing = 0.5 in Eq. [T Also the velocity law of a classi- 
cal ICastor. Abbott & Klein (1 19751) li ne-driven wind model can 
be approximated with /3 = 0.5 (c.f.. lGavlevll2000l) . Motivated 
by this, observational studies widely use the Eq. [T] and usu- 
ally consider /3 as a free parameter, yielding a typic al value of 
/3 = 0.7 - 1 for O star winds (e.g.. |Puls et al.ll 19961) . However, 



for some stars (e.g. , early B supergiants' ), higher /3 up to 3 may 
be determined (e.g. lCrowther et 31.11200 6*). 

Several variants of the simple formula Eq. [T] can be found 
in the literature and gene rally ai m to improve the fits of wind 
line profiles. For example. iHillien (j 1988.) introduced the velocity 
law that also encompasses the wind part clo s e to th e hydrostatic 
photosphere. Subsequentlv.' Hilher & Miller! (1 199 9*) introduced a 
second component in the beta law, which provides better fits to 
the blue edges of P Cygni fine profiles. This was also supported 
bv iGrafener & HamannI (l2005h . who found that different values 
of [3 describe the wind velocity law in different regions. The fail- 
ure of the simple /3-velocity law to describe the wind velocity 
law was also mentioned by Gayley et al. ( 1995). 

Although the original beta velocity law was derived theo- 
retically, there is no reason to believe that it provides the best 
approximation to the real situation. It is well known from hydro- 
dynamical simulations (IFeldmeier et al.l 1 19971 lOwocki & Pulsl 
11999 ) that the wind velocity field is variable with time, and con- 
sequently the beta law may at most describe mean flow. But even 
in stationary models, the wind velocity law may deviate from 
Eq. [T] at least because the wi nd is driven by different ions c lose 
to the star and at large radii (iVink et al.ll200UlKrtigkall2006l) . 

It would be highly desirable to develop formulae that are able 
to describe the velocity law in a more accurate form than Eq. [T] 
It is not clear how to develop a more general formula based on 
Eq. [T] Luckily, the problem of functional fitting is fully settled 
in mathematics, where the concept of orthogonal functions was 
introduced. Consequently, we propose to use the set of polyno- 
mials that may form an orthogonal basis in a given functional 
space, e.g.. 
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from which the terminal velocity follows as Voo — 51]"= i ^i- 
Equation |2a| can be equivalently rewritten in terms of Legendre 
polynomials. This guarantees that the derived coefficients de- 
crease with increasing order i and that the coefficients do no t 
significantly change with increasing n (e.g.. lMikulaieai2007h . 
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Fig. 1. Comparison of the fit to the calculated hydrodynamical 
model radial velocity of HD 249 1 2 (crosses) by the beta velocity 
law (solid line), and by the polynomial fit Eq.|2b]for n — 3. 



In this case instead of Eq. |2al the parameterizations is 



ViPi{x), with a; = 1 - — , 



(2b) 



1=0 



and the terminal velocity —J2"^iVi-To ensure the orthogo- 
nality in the interval [0, 1] w e use shifted Legendre po lynomials 
here. The first four are (e.g^ lKoornwinder et al.ll2010b 



Pa = 1, 



6x 



Pi = 2x - 1, 
1, and Pg = 20x^ - 30x^ + 12x - 1. (3) 



The relations between coefficients iii and Vi follow from Eqs.|2] 
and|3] 



Vo = -Vi + -V2 + -V3, 



V2 = -^V2 + -V3, 

6 4 
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Only three of these (vi) coefficients are independent. 

To demonstrate the convenience of the formulas Eq. |2] we 
fit the radial velocity structure obtained from hydrodynamical 
NLTE wind models with CMF line force of iKrticka & KubatI 
( 1201 Ol see Fig. [1). The best-fit parameters of the beta velocity 
law and polynomial fit Eq. |2b] with n = 3 for individual stars 
are given in Table [T] To test the global accuracy of individual 
orders of the velocity approximation Eq. |2] we plot the relative 
difference between the velocity from the numerical models 
and its approximations u^pp (Eq. |2b]for n — 1,2,3 and Eq. [U 
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averaged over individual stars in Table [T] 
as a function of radius in Fig.|2] 

Already the first term in Eq. |2al i.e. the approximation 
v{r) — vi{l — R^:/r) (in fact Eq. [T]for /3 = 1), provides rea- 
sonable approximation to the wind velocity law suitable for sim- 
plified analysis (see Fig.|2]i. Including the second term in Eq.l2al 
we obtain better approximation than the usual beta velocity law 
both for the velocity and its derivative using the same number of 
free parameters (/3 and Voo.p for the beta velocity law, vi and 
V2 for Eq. |2all. The agreement in the outer parts of the wind 
(t/ R* ^ 11) can be improved further by adding the third term 
in Eq.|2a|(see Fig.|2]i. 



Fig. 2. The mean relative difference (see text) between the ve- 
locity law (upper panel) and its derivative v' = dv/dr (bottom 
panel) taken from hydrodynamical models and its approxima- 
tions averaged over stars from Table[T]as a function of radius. 



It has not escaped our attention that even relatively high or- 
der polynomials do not provide reasonable fits in the region close 
to the star, for r/i?, < 1.1. We do not aim to fit the velocity in 
this region, because here the wind density approaches the hy- 
drostatic density stratification. Consequently, it would be better 
to derive the velocity from the continuity equation using density 
from the hydrostatic equilibrium equation instead of a polyno- 
mial. 

Many studies fit the observed line profiles of particular 
stars usi ng sophisticate d radiative tra nsfer calculations (e.g., 
iPuls et al . 1996, Hillier & Milleilll999l) by assuming the veloc- 
ity law of type Eq. [1] and trying to find the most suitable value 
of the parameter f3 in combination with the value of the termi- 
nal velocity Voo- Although our formulae Eq.|2]were derived to fit 
the results of numerical hydrodynamical calculations, it might be 
interesting to use them to also fit the beta velocity profiles. The 
test showed that the polynomial approximation Eq.|2]with n = 2 
provides very good fits to beta laws (Eq. [TJ bot h with small < 
land large /3 > 1 and to the combined law of iHillier & Millerl 
(Il999h . Only for large /3 > 2.3 is it better to increase the poly- 
nomial degree (n = 3) to assure the positivity of the fit close to 

the star The modified beta law v(r) = Vqc,^ (l ~ ^^r)^ with 
additional parameter 6 < 1 ( Puis & Hummed 1988i) is frequently 
used to improve the fit close to the star. It can be satisfactorily re- 



produced by either introducing a fixed low value of vq in Eq.l2al 
or by a substitution x = 1 — bR^,/r in Eq.jjb] Moreover, our ex- 
periments with fitting the trial P Cygni line profiles have shown 
that the fits using Eq.|2]and the beta velocity law require a com- 
parable number of line profile evaluations. A reasonable fit strat- 
egy could be to derive Vao from observations and then keep the 
sum -Li = ^oo fixed. 

We conclude that the beta velocity law provides reasonable 
approximation to the hydrodynamical wind velocity law and its 
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Table 1. Best-fit parameters to the radial velocity from hydrodynamic models (iKrticka & KubatlbOlOl) for selected stars and indi- 
vidual velocity parameterizations. 



Star 


HD number 








^lawEq.UJ 




Eq.|2bjfor n = 3 








R, 


M 


Tcff 




P 


vo 


Vl 










[Rq] 


[Mq] 
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O stars 


^Per 


HD 24912 


14.0 


36 


35 000 


2130 


0.70 


1270 


980 


-220 


60 


lOri 


HD 37043 


21.6 


41 


31400 


2170 


0.83 


1190 


1040 


-120 


30 


15 Men 


HD 47839 


9.9 


32 


37 500 


3240 


0.90 


1700 


1590 


-90 


20 




HD 54662 


11.9 


38 


38 600 


2200 


0.79 


1230 


1070 


-180 


-10 




HD 93204 


11.9 


41 


40 000 


2210 


0.79 


1240 


1070 


-170 





COph 


HD 149757 


8.9 


21 


32 000 


1940 


0.85 


1060 


930 


-50 


70 


68 Cyg 


HD 203064 


15.7 


38 


34 500 


1970 


0.76 


1120 


920 


-150 


50 


19Cep 


HD 209975 


22.9 


47 


32 000 


2090 


0.82 


1160 


1010 


-120 


30 


Central stars of planetary nebulae 


NGC 2392 




1.5 


0.41 


40000 


540 


0.86 


290 


270 


-50 


-20 


NGC 3242 




0.3 


0.53 


75 000 


2070 


0.68 


1260 


920 


-290 


60 


IC 4637 




0.8 


0.87 


55 000 


1350 


0.89 


710 


680 


-70 


-30 


IC 4593 




2.2 


1.11 


40000 


730 


0.88 


390 


360 


-60 


-30 


IC418 




2.7 


1.33 


39 000 


800 


0.87 


420 


380 


-70 


-20 


Tc 1 




3.0 


1.37 


35 000 


910 


0.95 


460 


440 


-30 


-10 


NGC 6826 




2.2 


1.40 


44000 


880 


0.88 


470 


430 


-60 


-30 



Notes. The stellar parameters (radius mass M, and the effective temperature Teff) are taken from lRepolust et aT] ( l2004h . lMarkova et alj ( l2004h . 
[Pauldrach et al.. (.2004.') . and .Martins et aL(.2005,) . 

derivative. The approximation can be improved with polynomial 
fits Eq. |2a] (or, equivalently Eq. l2b] i using the same number of 
free parameters as in the commonly used beta velocity law. 
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